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abstract 



The present thesis will allow Users to analyze rectangular 
plates with or without a rectangular hole with its edges parallel 
to the edges of the plate , by using an approximate grid analogy. 

The method is completely computerized and the user need 
only specify properties, geometry, loading and kind of output de- 
sirable. The computer programs then will generate a deck of cards, 
that form a complete STRUDL input. STRUDL itself is an equivalent 
of STRESS but with much more useful capabilities. 

Initially it was intended to make everything an one step 
process, i.e. to read input in a free-field form, generate all 
necessaxy statements of STRUDL, store them on a file and let 
STRUDL read the file and give the desired output. However due to 
several delays in the Civil Engineering system Laboratory, STRUDL 
is not yet available and CDL (Command definition language) is now 
in the process of experimentation. In consequence the user has to 
supply his input according to specific format and the linking be- 
tween this subsystem and STRUDL is not automatic. In other words 
the process is a two-step one and the user has to take the output 
of this program and use it as a STRUDL job. 

The solution to the flexural plate problem contained herein, 
is extremely general regarding boundary conditions and loadings 
so that almost eveiy engineering problem of this nature can be 



taken care of 



The particular connection of tho present job with Naval 
Architecture is that solutions to hotter, and side plating with 
hydrostatic loads* bulkheads with openings and decks with cargo 
holes will become readily available with a minimum of human effort 



and intervention* 
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CHAPTER I 

A ?iT3CDUCglO» AKD JCOFE 

The development of tine large digital computers was a deci- 
sive step towards the solution of many engineering problems. Indeed 
during the past few years people looked back and rcoxa 

but it was impossible to carry it out by hand, xx, was a 
simple matter to solve all these problems then by a digital com- 
puter. 

Among all those other problems one is the structural analy- 
sis of plates a special case of which is the rectangular plate in 
flexure. Several scientists and engineers have done a marvellous 
job In solving special casos of it. 

It seems that there would be two distinct approaches in 
coding the solution of the problem into a computer: 

(i) Take each special case and code its exact solution. 

(ii) Invent a unified approach applicable in all or most 
special cases. 

The first approach requires a tremendous aremount of work. 
Besides it could never be called general. For, no matter how many 
special cases one solves, there will always be some other special 
case left out. 

The second approach seems quite attractive, with the only 
disadvantage that it is an approximate one. However, in casos where 
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an exact solution is not known, an approximate one is best than 
nothing. Besides, very of ton engineers are interested in getting 
a first estimate of what is happening in a system. 

In the structural analysis of rectangular platos, the grid 
analogy which will be deccribod in detail in the next chapter, is 
among the best approximate methods to attack tho general problem. 

In addition since several methods to solve skeletal structures 
have now boon developpod, ono has only to transform the plate 
into such a structure, that existing computer formulations can 
solve. 

The goal set forward by the author has been to develop 
a set of computer programs, able to reduce ary rectangular plate 
to its skeletal equivalent and then solve it as such. Although 
the sample problem included in the present work is a very simple 
one, the capabilities of the programs developed aro quite general. 
In fact the choice of the sample problem, made strictly because 
its solution already exists, is a very unfortunate one in illustra- 
ting the application of the method into the more general problem. 

The case of a rectangular hole on the plate was included 
here, due to tho interest of Naval Architects on the effect of 
cargo openings on the strength of ship's decks. Also man-holes 
on vertical bulkheads forming boundaries of fu oil -oil tanks con- 
stitute another field of application. 

What really makes the method powerful is its generality 
regarding the boundary conditions and the loading. 
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Indeed, point supports at any point over the area of the 
plate, line boundary conditions over all or any part of the edges, 
concentrated point loads, uniform or linearly distributed line 
and area loads or any combination of them, are all cases that can 
be taken car© by the programs with a minimum of effort. In addition 
the ability to get ary kind of output desired might be very useful 
in some cases. 

In the following chapters the analogy will bo described 
and the method of its computer application will be explained. In 
addition directions will be given to the usors on how to use the 



programs. 
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CHAi-TEa II 

jmA m& x. 



2*1*- general 

Perhaps one could say with complete Justification that the gene- 
ral problem of a plate In flexure* has been solved* at least in principle* 
Indeed the differential equation describing the phenononon is known for 
many years and it is only necessary to consider the loading and boundary 
conditions applicable in each particular problem. However this last step, 
namely* the solution of the differential equation for a given set of 
boundary conditions and a given loading, eonetines presents innumerable 
difficulties from the analytical point of view* Sven in the case whore 
a solution can be found in a reasonable number of steps, it is very 
often expressed in the form of an infinite series, which especially 
in the case of moments and ah *ars converge rather slowly* a fact that 
reduces thoir value. 

It is the above difficulty that made engineers search for other, 
perhaps less accurate but more easily developed solutions of the plate 
problem* There have been three approaches that yield approximate solu- 
tions to the plate stresses problem. 

(i) Finite differences methods. 

(ii) Grid analogies. 

(iii) Finite element methods. 

Cf these, the finite differences approach, depending upon the 
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boundary conditions, very ofton requires fictitious extensions of 
members or consideration of members with fictitious properties. This ma- 
kes a rather complicated task the coding of the general rectangular plate 
solution to a digital computer. There is also some suspicion that this 
method gives results deviating from the exact ones (those that they 
exist) more than those computed by the other two methods. 

The finite element method seems to be quite promising not only for 
tho rectangular plate case but for shells of arbitrary shape. Conside- 
rable numbor of people is working on that method, in the civil Engineering 
department of H.I.T* , bat there will be some time before this method be- 
comes available to users. 

The grid analogy was selected by the author for tho following 

reasons. 

(i) It is rather simple and gives very satisfactory results 
regarding moment intensities and deflections as it was indicated in (1), 
(2) and (5). 

(ii) It can be rather easily programmed to a digital computer. 

(iii) The solution to the equivalent grid can be very easily- 
found by the existing STRESS program or by its new equivalent STBUDL, 

( 7 ), ( 8 ), ( 9 ). 

(iv) The present work intended primarily for the field of Naval 
Architecture as well as other professions interested in structural 
engineering deals almost exclusively with rectangular flat plates or 
plates with small curvature which can be neglected, without considerable 
reduction in accuracy. 
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2*2. Description of grid analogy 

The grid analog of a flat plate has been represeftted differently 

by different people. Lightfoot in (3) uses a square gridwork consisting 
of grid beams with flexural rigidity and a set of torsion bars with tor- 
sional rigidity with no diagonals. Efimba in (6) in a program written 
to analyse bridge docks only, uses TambergSs method of analysing beam 
and slab type bridges as torsicnally weak grids (10). 

The approach taken by Erennikoff in (2) and Yettram and Husain in 
(1) is, that the plate can be represented by a plane structure composed 
of articulated bars, joined by pins normal to the plane of the structure. 
The pattern has to repeat itself throughout the area of the plate and 
the outline and external dimensions of the plate should be the same 
with the corresponding ones of the equivalent structure. Repetition of 
the pattern does not however mean that the size of the unit has to re- 
main the same. In fact by changing the size of the unit we can satisfy 
any existingdeviations from the simplest case of a simple rectangular 
plate, like differences in thickness, holos or differences in material 
properties. 

Erennikoff sot up the problem in a very detailed form and 
developed the equations that will define the properties of the elements 
of the unit. The restrictions that he imposes in order to solve for 
the properties of the members are, that the frame work has to deflect 
properly: 

(i) When subjected to pure bending in the X direction with 
no bending deformations in the Y direction. 

(ii) When in pure bonding in the Y direction with no deforma- 
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tion in the X direction. 

(iii) When subjected to pure torsion on the X and Y planes. 

These are general restrictions applicable in any case and are 
useful in examining the suitability of any particular unit to be 
investigated. The third condition falls out whenever the pattern has 
the same form in the direction of both axes. 

Subsequently Hrenikoff investigated two patterns as shown in 
fig's (2.1) and (2.2) and proceeded to develop their properties. 

Of these two patterns the first one restricts the solution 
to values of Poisson's ratio of 1/3 only. 

The second one although valid for 
any Poisson's ratio has the disadvan- 
tage that is a square one. This causes 
a serious problem in the case that 
we cannot satisfy the boundary don- 
ditions with a square pattern. We 
also notice that for each square we 
have 10 slender members, which would 
make a hand-solution more difficult 
or a computer solution longer. 

A rectangular pattern with diagonals 
only is much more general and simple. 
This type of pattern has been very 
recently investigated by Yettram 
and Husain in (1). Their pattern 
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looks to be the boot of the existing ones because it is simple, easily 
developed and considerably general. Due to the fact that the net hod 

is relatively new and capes iaily useful in other kinds of problems 
(those satisfying the Langrangian equation like thcmal stresses in 
a plats ox special cases of fluid flow), it is believed that a rela- 
tively detailod discussion ot it must be given. Ibis will be doae in 
the next chapter. 

2.3. 

The idea of the aodaL is that w* can substltudc a rectangular 
network consisting of slendor members only, in place of a particular 
plate element. Fig. (2.3a) shows a rectangular element of plating and 
Fig* (2.3b) its equivalent grid element. The unknown properties of the 
grid element members are the flexural rigidity of the side beans I c 
and I s , the torsional rigidity of the side beams — p and G ^ c 

and the flexural rigidity of the diagonals . That gives 5 unknown 
properties i.e. I c » 1 8 * Id , — and — . 

The equations giving the above 5 properties can be developed 
by equating the rotations at the nodes of the grid element with those 
of the plating element of equal size. Although the method doasn*t give 
the cross-sectional area of the beams, it gives all the necessary in- 
formation to make tho solution of the grid frame-work possible. Of 
course tho plate ard the grid elements nuat bo subject to statically 
equivalent torques arid moments. 

Considering first the plating element we examine three cases. 

(i) Equal and opposite bonding moment intensities on the x sides 
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(ii) Equal and opposite intensities on y sides . 

(iii) Self-equilibrating torque intensities along 4 edges. 

Fig's (3d)* (3«) and (3g) show clearly the moment and torque 

distributions and the corresponding angles. The values of the angles 
in terms of the moment intensities can be found by differentiating the 
expressions for the deflections, (ll). This was done In (1) where 
they came up with the following expressions for the angles. 



0i - 






( 2 . 1 ) 



- 



EC W 1 / > 1) 



e 



H \A ^ Mz 



E (V*/'*) 



( 2 . 2 ) 

(2.3) 

(2.4) 



^ 3 h ( 1 + r-) 
e Cw 3 In) 



(2.5) 



For the grid element loaded a3 shown in Fig’s (3d), (3f) and 
(3h). Ref (1) gives the flexibility matrix (shown here in Fig 2. 4) 
and concludes that only 2 out of the 12 equations are independent when 
wo have only side bending moments or self-equilibrating torques. 

By selecting properly these 2 equations in each case we come up with 
the following expressions for the grid angles. 
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Ity equating the rotations for the plate and the grid elements 
wo can get the following equations giving the required properties of 
the grid member. [The method is explained in nor© detail in (1)J . 
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Where 



k is the ratio of the side beans 

|X is the Poisson’s ratio 

h is the thickness of {.dating 

o is the loung’o modulus of elasticity 

G is the modulus of rigidity 

r is tii# diagonal to side boa.; ratio 

For future use we will transform the above relations by using* 
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into toe following fora* 
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le a constant for a given isotropic plate of constant thickness. 

The potentialities of the above equations aro Ian od lately 
apparent. They can be used for the solution of ary plating problem, 
with ary boundary conditions and load and any number of holes provided 
that we can form a grid fine enough so that it will represent a parti- 
cular plating element, i.e. provided that we will satisfy equations 
( 2 . 16 ) through ( 2 , 23 ) at each particular element in the division of 
the plate. Also by forming a proper grid we can take caro of non- 
isotropic plates. 



